Transport of fullerene molecules along graphene nanoribbons 
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We study the motion of Cgo fullerene molecules (buckyballs) and short-length carbon nanotubes on graphene 
nanoribbons. We demonstrate that the nanoribbon edge creates an effective potential that keeps the carbon struc- 
tures on the surface. We reveal that the character of the motion of Ceo molecules depends on temperature: for 
low temperatures (T < 150 K) the main type of motion is sliding along the surface, but for higher temperatures 
the sliding is replaced by rocking and rolling. Modeling of the buckyball with an included metal ion, such as 
Fe~ @C60, demonstrates that this molecular complex undergoes a rolling motion along the nanoribbon with the 
constant velocity under the action of a constant electric field. The similar effect is observed in the presence of 
the heat gradient applied to the nanoribbon, but mobility of carbon structures in this case depends largely on 
their size and symmetry, such that larger and more asymmetric structures demonstrate much lower mobility. Our 
results suggest that both electorphoresis and thermophore sis can be employed to control the motion of carbon 
molecules and fullerenes and, for example, sort them by their size, shape, and possible inclusions. 



I. INTRODUCTION 

Over the past 20 years nanotechnology has made an im- 
pressive impact on the development of many fields of physics, 
chemistry, medicine, and nanoscale engineering [1]. After the 
discovery of graphene as a novel material for nanotechnol- 
ogy [2], many properties of this two-dimensional object have 
been studied both theoretically and experimentally [3, 4]. In 
addition to its unusual properties, graphene may play an im- 
portant role as a component of more complex systems of nan- 
otechnology. In particular, atoms or molecules deposited on a 
surface of single-layer graphene [5] are expected to move with 
just two degrees of freedom, demonstrating various phenom- 
ena that can be employed for nanomechanics and nanoscale 
molecular motors [6]. 

Numerical modeling of the motion of atoms and molecules 
on a graphene sheet was considered very recently in a num- 
ber of studies [7-9]. In the first approximation, one may 
treat graphene as a planar two-dimensional substrate that cre- 
ates an effective two-dimensional periodic potential for ab- 
sorbed molecules. Stochastic motion of noble gases on such 
a substrate potential was analyzed in Ref. [7]. It was shown 
that Xe atom is trapped in a potential well of the substrate 
at high friction while He atom can freely diffuse. Diffusive 
motion of Cgo molecules on a graphene sheet was modeled 
in Ref. [8]. Analytical and numerical studies demonstrated 
that the translation motion of Ceo molecule near a graphene 
sheet is diffusive in the lateral direction, whereas in the per- 
pendicular direction this motion can be described as diffusion 
in an effective harmonic potential (bounded diffusion). It was 
also shown that the motion of Ceo over graphene sheet is not 
rolling. More detailed analysis of the structure of the effec- 
tive two-dimensional periodic potential for Cgo molecules was 
carried out in Ref. [9]. It was shown that the energetically fa- 
vorable motion of fullerene molecules is sliding along the sur- 
face when one side of the Ceo molecule remains parallel to the 
plane of the graphene and the molecule center is placed above 
the valent bond of the nanoribbon. However, all those results 
were obtained for the case when the graphene nanoribbon is 



assumed to be flat and rigid, so that the motion of its atoms is 
not taken into account. 

Directed motion of Ceo molecules along the graphene sheet 
under the action of a heat flow was studied numerically in 
Ref. [10]. In their calculations, the authors employed the so- 
called Nose-Hoover thermostat, which however is known to 
produce incorrect results in the nonequilibrium molecular dy- 
namics simulations [11, 12]. In addition, the directed motion 
of the fullerenes were not analyzed. 

By now, novel nanocomposite structures consisting of a 
stacked single graphene sheet and Ceo monolayer were pro- 
duced experimentally [13]. It was found that the intercalated 
Ceo molecules can rotate in between graphene sheets. Such 
a rotation of a single Cgo molecule on a sheet of graphene 
should result into its rocking and rolling motion, in contrast 
to the statement of Ref. [8]. In order to resolve those contro- 
versies and model more general dynamics closely resembling 
the experimental conditions, one have to take into account the 
motion of carbon atoms of the nanoribbon that will modify 
its planar geometry due to the interaction with Cgo fullerene 
molecules. 

In this paper, we study numerically different types of mo- 
tion of carbon structures, such as Cgo fullerene molecules 
and short-length carbon nanotubes, on finite-width graphene 
nanoribbons. We employ molecular dynamics simulations 
taking into account the motion of all carbon atoms, and 
demonstrate that an edge of the nanoribbon creates a poten- 
tial barrier for the absorbed nanostructures supporting their 
transport along the nanoribbon. We study different types of 
dynamics including stochastic Brownian like motion as well 
as directed motion of fullerenes driven by the electric field 
and temperature gradient. We reveal that the motion of Ceo 
fullerene molecules depends on the temperature, so that for 
low temperatures the main type of motion is sliding along the 
surface, but for higher temperatures the sliding is replaced by 
rocking and rolling. For the directed motion, the mobility de- 
pends largely on the shape and symmetry of the carbon struc- 
tures, so that larger and asymmetric structures demonstrate 
much lower mobility. Our results suggest that both elector- 
phoresis and thermophoresis can be employed to control and 
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FIG. 1 : Motion of (a) fullerene Ceo molecule, and (b-d) shot-length 
(5,5) carbon nanotubes Cioo, (c) Ciso, and (d) C260, respectively, 
along a finite-width zigzag carbon nanoribbon. Valent bonds of the 
nanoribbon are shown in green, while the valent bonds of the ab- 
sorbed fullerene-like molecules, are shown in red. Temperature is 
T = 300K. 



direct carbon structures and fullerenes. 



II. MODEL 

The first study of the van der Waals (VDW) energy for 
the interaction of fullerene molecule Ceo and substrate with 
atomic surface corrugation was presented in Ref. [14]. The 
VDW energy was calculated using linear response theory to 
evaluate the dipole-dipole interactions between the molecule 
and the substrate. Interaction of the fullerene Ceo with 
graphite surface was described by the 12-6 Lennard- Jones 
(LJ) potential 

V(r) = C 12 /r 12 - C 6 /r 6 (1) 

with the coefficients C 12 = 12000 eVA 12 and C 6 = 
15.2 eV A 6 . 

Experimental evaluation of the interaction of molecule C60 
with single-wall carbon nanotubes (SWNTs) and graphite was 
obtained in Ref. [15], on the basis of the thermal desorption 
spectroscopy. It was shown that these interactions can be 
modeled, with a high accuracy, by the LJ potential (1) with 
coefficients C 12 = 22500 eVA 12 and C 6 = 15.4 eVA 6 . This 



potential has been employed for the modeling absorption of 
Cgo molecules by single-walled carbon nanotubes with open 
ends. 

Detailed modeling of the motion of a Cqo molecule inside 
a single-walled carbon nanotube was carried out in a series 
of studies [16-19]. In particular, Ref. [16] studied the mo- 
tion of C 60 inside (11,11) SWNT (the system C 60 @SWNT), 
Ref. [17] analyzed the oscillatory motion between the edges 
of a finite-length (10,10) SWNT, and Ref. [18] demonstrated 
that the energy dissipation of the oscillatory motion of C^o de- 
pends substantially on the radius of SWNT and the presence 
of defects. Thermophoretic motion of buckyball enclosed in 
a series SWNTs was analyzed in Ref. [19]. It was shown that 
inside the nanotube the encapsulated fullerenes move in the 
direction opposing the thermal gradient, i.e. in the direction 
of the thermal flow. Similar behavior of the Cqq molecule was 
predicted for a graphene sheet where the fullerenes was shown 
to move in the direction of a heat flow [10]. 

In this paper, we study the dynamics of a molecular sys- 
tem consisting of a fullerene molecule Cgo or short-length car- 
bon nanotubes placed on a graphene nanoribbon, as shown in 
Figs, l(a-d). We employ numerical molecular-dynamics sim- 
ulations taking into account the motion of all carbon atoms. 
Inside of the fullerene molecule and a sheet of graphene the 
molecules two neighboring carbon atoms create a valent bond 
[see Fig. 2 (a)]. We describe the energy associated with the 
deformation of the valent bond created by two carbon atoms 
by the following interaction potential 

U 1 {p)=e 1 {eM-a(p-pi)]-l}\ (2) 

where p is the length of the valent bond, t\ = 4.9632 eV and 
pi = 1.418 A is the energy and the equilibrium length of 
the bond, parameter a = 1.7889 A -1 . Each carbon atom is 
placed at the vertex of three planar valent angles. The corre- 
sponding deformation energy of a plane valent angle created 
by three atoms [see Fig. 2 (b)] can be described by the inter- 
action potential of the form, 

U 2 (<p) -e 2 (cos^ + l/2) 2 , (3) 

where ip is value of the valent angle (equilibrium value of the 
angle ip = 2n/3), energy e 2 = 1.3143 eV. Each valent bond 
is simultaneously belong to two planes, so that the deforma- 
tion energy of the angle formed by two such planes can be 
described by the following interaction potential, 

U i (<j>)=e i (l-cosct>), (4) 

where - is the corresponding angle (in equilibrium, <j> = 0), 
and the index i = 3, 4, 5 describe the type of the dihedral 
angle - see Fig. 2 (c), (d), (e). Energy £3 = £4 = 0.499 eV, 
£5 <C £4, so that the latter contribution to the total energy can 
be neglected. 

More detailed discussion and motivation of our choice 
of the interaction potentials (2), (3), (4) can be found in 
Ref. [20]. Such potentials have been employed for modeling 
of thermal conductivity of carbon nanotubes [21, 22] graphene 
nanoribbons [20] and also in the analysis of their oscillatory 
modes [23-25]. 




FIG. 2: Typical configurations of the nanoribbon bonds containing up to i-th nearest-neighbor interactions for: (a) i = 1, (b) i = 2, (c) i — 3, 

(d) i = 4, and (e) i = 5. 



We describe the nonvalent interaction of carbon atoms 
belonging to the fullerene molecule and graphene by the 
Lennard- Jones 12-6 potential 



U LJ (r) = 46 
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where r - distance between the atom centers, interaction en- 
ergy Co = ecc = 0.002635 eV, and parameter a = tree = 
3.3686 A. Potential (5) coincides with the potential (1) for 
the coefficients C 12 = ^cc^cc = 22500 eV/A 12 , C 6 = 
4e cc cr6 c = 15.4 eV/A 6 (see [15]). 



III. APPROXIMATION OF RIGID MOLECULES 

First, we study the energy of the coupled system "Ceo + 
nanoribbon" as a function of the relative position of the buck- 
yball on the surface. We consider a square sheet of graphene 
of the size 4.1 x 4.1 nm 2 [see Fig. 3 (a)] and assume that both 
the Ceo molecule and the graphene are in the ground state and 
all valent bonds are not stretched, and the angles between the 
bonds do not deform. Then the energy of the complex is de- 
fined by a sum of all nonvalent interactions (5) between the 
carbon atoms of the fullerene molecule and the graphene sub- 
strate. 

Without restrictions of generality, we assume that the sheet 
of graphene remains at rest in the (x, y) plane, so that the po- 
sition of the Ceo molecule is defined by the coordinate (x, y) 
of its center on the plane and the distance z from the center to 
the substrate surface. We fix the center coordinate (x, y) and 
minimize the energy E of the whole molecular complex as a 
function of the distance z and all angles that define the ori- 
entation of the Ceo molecule in space. Then we obtain two- 
dimensional function E(x, y) that describes a change of the 
energy of the system when the fullerene molecule changes its 
position on the graphene sheet. 

The effective potential energy E(x, y) shows that in the 
ground state the Ceo molecule takes the position above the 
center atom of the graphene sheet at the distance zq = 
0.635 nm from the surface, so that the hexagonal face closest 
to the surface is parallel to it. Figure 3 (c) shows the depen- 
dence of the energy E on the distance z to the surface. As 
follows from this dependence, in order to escape the surface 
the molecule requires the energy AE = 0.73 eV. In the ear- 
lier study [8], similar values were defined as zq = 0.65 nm 
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FIG. 3: (a) Energy E of the fullerene Ceo molecule as a function of 
its position on a square sheet of a single-layer graphene of the size 
4.1 x 4.1 nm 2 . Red lines mark the valent bonds of the graphene, 
where the point [x, y) defines the projection of the molecule center 
on the plane, (b) Shape of the potential energy E(x, y) at the center 
of the graphene sheet, circles mark the position of the carbon atoms 
of the graphene. (c) Interaction energy as a function of the distance 
z between the center of Ceo molecule and the plane of graphene. 
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FIG. 4: Change the energy of Ceo molecule due its rotation around 
the edge closest to the graphene sheet. Curve 1 : rotation from one 
hexagon face to the neighboring hexagon face; curve 2: rotation from 
the hexagon face to the neighboring pentagon face. 



and AEq = 0.76 eV, and the experimental evaluation of the 
coupling energy [15] gave the value AEo = 0.85 eV. 

Dependence of the energy E of the coupled molecular sys- 
tems on the position coordinate (a;, y) of the molecule center 
on the plane of the graphene sheet is shown in Figs. 3 (a,b). 
At the center of the sheet, the energy depends on the position 
of the fullerene molecule, and the positions above the carbon 
atoms of the substrate are more favorable and they correspond 
to local minima. The energy profile suggests that the molecule 
may slide along the graphene sheet with its hexagonal face 
parallel to the surface following the structure defined by the 
valent bonds of the substrate. This motion would require over- 
taking the energy barrier AEi = 0.0015 eV, whereas the 
sliding motion in other directions would require higher en- 
ergy to overtake the barrier AE 2 = 0.015 eV. We notice that 
the sliding of the Cgo molecule between two parallel sheets 
of graphene was studied earlier in Ref. [9], where was found 
that the corresponding barriers are: AE\ = 0.002 eV and 
AE 2 = 0.02 eV. 

When the molecule approaches the edge of the nanoribbon, 
the total energy of the system E(x, y) grows, as shown in 
Fig. 3 (a), and at the edge of the nanoribbon potential barrier 
AE% > 0.3 eV repels the molecule from the edge. 

In addition to the sliding, the Cgo molecule may also roll 
along the surface. To evaluate when this rolling becomes pos- 
sible, we analyze a change of the energy of the ground state 
when the molecule rotates around the edge of its hexagonal 
face closest to the graphene surface. 

In the ground state the hexagonal face of the fullerene clos- 
est to the graphene sheet is parallel to the substrate surface. If 
we fix one edge of this face and rotate the buckyball around 
it modeling the rocking of the Cqq molecule, we can find the 
dependence of the energy E on the rotation angle <fr, shown in 
Fig. 4. For such a procedure, two cases are possible. In the 
first case, the molecule rotates replacing one hexagonal face 
by another similar face, the energy E(<p) for this rotation has 
the form of a symmetric two-well potential, as shown in Fig. 4 
(curve 1). For the other type of rotation, the hexagonal face 




FIG. 5: Stochastic dynamics of the fullerene molecule Cbo placed on 
a graphene sheet with the size of 4.1 x 4.1 nm 2 . (a,c) Trajectory of 
the molecule center shown as a projection on the plane of graphene 
at the temperatures T = 30K and 300K, respectively. Red color 
shows the valent bonds of carbon atoms of the graphene, red dots 
show the fixed atoms. (b,d) Temporal evolution of the z projection 
of the center of the hexagonal face relative the molecule's center of 
mass, for T = 30K and 300K, respectively. 



is replaced by a pentagonal face also parallel to the substrate 
surface, in this later case the energy has the form of asymmet- 
ric two- well potential also shown in Fig. 4 (curve 2). In both 
these cases of the molecule rolling, the molecule should over- 
take the energy barrier AE± = 0.052 eV substantially higher 
than the energy barrier AE 2 required for its sliding. 

Our analysis suggests that the Ceo molecule may create 
long-lived coupled states with the graphene nanoribbons or 
graphene sheets of a finite extent, and its sliding along the sur- 
face is more energetically favorable than its rocking or rolling. 
Thus, we may expect that for low temperatures the main type 
of motion of the fullerene molecules would be sliding, and it 
will be replaced by rolling at higher temperatures. To verify 
these ideas, below we provide the numerical modeling of the 
motion of the fullerene molecules along the graphene surface 
taking into account all degrees of freedom of the coupled sys- 
tem and considering different values of temperature. 



IV. STOCHASTIC MOTION OF Cm ON A GRAPHENE 
SHEET 

First, we study the stochastic motion of the fullerene 
molecule Cqo placed on a thermalize sheet of graphene. We 
consider a flat square sheet of graphene of the size 4.1 x 
4.1 nm 2 composed of 678 carbon atoms in the (a;, y) plane, 
as shown in Figs. 5(a,c). To keep the graphene in the plane, 
we fix the four edge atoms, shown by solid red circles in the 
corners of the plots of Fig. (a,c). We place the molecule Cgo 
at the middle of the sheet and place the whole system "buck- 
yball + graphene" consisting of N = 678 + 60 atoms into 
Langevin thermostat modeled by stochastic forces. The mo- 
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tion equations are taken in the form 

3H 

M l ix l = -- TMiUi + Ei, (6) 

dui 

where H is Hamiltonian of the molecular complex, the indices 
i = 1, 2, N stand for the atom number, the vector = 
(v,i i,Ui2, u i,3) describes the position of the atom, and Mj is 
the atom mass. Langevin collision frequency is T = l/t r , 
where t r - particle relaxation time, and Ej = (& i,£»,2j £.1.3) 
is a three-dimensional vector of Gaussian stochastic forces de- 
scribing the interaction of the i-th atom with thermostat. 

We consider a hydrogen-terminated graphene, where all 
edge atoms correspond to the CH group. In our model, we 
treat such a group as a single effective particle at the loca- 
tion of the carbon atom, so that we take the mass of atoms 
inside the graphene sheet and fullerene atoms as Mj = 12m p , 
but for the atoms at the graphene edge we consider an effec- 
tive particle with a larger mass Mj = 13m p (where m p — 
1.6603 x 10 -27 kg is the proton mass). 

To analyze the thermalized dynamics of the system, we in- 
troduce a random white noise with the correlation functions 

<6,*(ti)&,i(t2)> = 2M i Tk B T5 ij S kl 5(t 2 - h), (7) 
i,j = l,2,...,N, k,l = 1,2,3, 

where T is the temperature of the Langevin thermostat. We 
select relatively large relaxation time t r = 1 ps to reduce the 
effect of viscosity in the molecular dynamics, and integrate 
the system of equations (6) for the time interval t = 4 ns (for 
the period t = I0t r = 10 ps the system should relax into the 
equilibrium state with the thermostat and then it will evolve 
as a fully thermalized system) 

Our numerical simulations demonstrate that for the tem- 
peratures T < 500K the molecular system preserves its ge- 
ometry, namely the buckyball Cgo remains attached to the 
graphene sheet, and it undergoes stochastic motion not ap- 
proaching the edge of the sheet. For very low temperatures 
T < 150K the main type of motion is sliding along the sur- 
face, while for higher temperatures the sliding is replaced by 
rolling. Figures 5(a,b) show an example of such a motion at 
T = 30K when the buckyball moves along the valent bonds 
of the graphene sheet [see Fig. 5 (a)] changing time to time 
its sliding face [see Fig. 5 (b)], such that for the time period 
of 4 ns it changed the face 3 times, as shown in Fig. 5(b). 
For higher temperatures such as T = 300K the dynamics is 
clearly stochastic, and the motion does not reflect much the 
discrete structure of the graphene substrate, with the contact- 
ing face being changed continuously [see Figs. 5(c,d)]. We 
observe that in this case the main type of motion is rolling 
along the graphene sheet, which is energetically more favor- 
able than sliding. 



V. ELECTORPHORESIS OF FE @ C 60 

In order to model numerically the motion of a spherical C@o 
molecule along the graphene nanoribbon under the action of a 




FIG. 6: Transport of the molecular complex Fe~@Ceo along the 
nanoribbon with the dimensions 157 x 2 nm 2 , under the action of the 
constant electric field E = 1CT 7 V/m, at temperature T = 300K. (a) 
Evolution of the coordinate of the molecule center for four different 
thermalizations of the nanoribbon (curves 1,2,3, and 4). An auxiliary 
straight line show the slope of the curves corresponding to the max- 
imum velocity of the molecule, v m — 390 m/s. Dashed straight line 
marks the location of the edge of the nanoribbon. (b) Temporal evo- 
lution of the z projection of the center of the hexagonal face relative 
the molecule's center of mass, corresponding to the curve 1 in (a). 



constant electric field, we include an ion of iron Fe into the 
buckyball, and consider a composite structure Fe~@C6o- 

We describe the interaction of the iron molecule with the 
carbon atoms of the molecule Cgo by employing the LJ poten- 
tial (5) with the interaction energy eo = £FcC = 0.0031 eV 
and a = rFec/2 1 ^ 6 , where the equilibrium distance between 
the atoms is rp G c = 3.51 A. We study the motion of the 
molecular complex Fe~@Ceo along the zigzag nanoribbon 
consisting of 640 segments with 20 carbon atoms in each. 
Such a ribbon has the dimensions 157 x 2 nm 2 and it consists 
of 12800 carbon atoms (a part of the nanoribbon is shown in 
Figs. 1 (a-d). 

We place an undeformed nanoribbon in the (x, y) plane 
along the axes x and fix it in the plane by pinning its four edge 
atoms. Then, we place the system into the Langevin thermo- 
stat by solving numerically the system of stochastic equations 
(6) with T = l/t r , where t r = 1 ps. During the time t = 10i r 
the system will get thermalized and then we add the molec- 
ular complex Fe~@C6o- By applying the constant electric 
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FIG. 7: Transport of the Cgo molecule along the nanoribbon with the 
size 157 x 2 nm 2 , under the action of a heat flow, (a) Evolution of 
the coordinate of the molecule center for four different thermaliza- 
tions of the nanoribbon (curves 1,2,3, and 4). An auxiliary straight 
line shows the slope of the curves corresponding to the maximum 
velocity of the molecule, v m — 190 m/s. (b) Temperature profile 
along the nanoribbon, a solid straight line corresponds to the tem- 
perature gradient of dT = 0.30 K/nm. Grey shaded areas show the 
edge areas of the nanoribbon subjected to two Langevin thermostats 
at T+ = 330 K and T_ = 270 K, respectively, (c) Temporal evo- 
lution of the 2 projection of the center of the hexagonal face relative 
the molecule's center of mass, corresponding to the curve 1 in (a). 



FIG. 8: (a) Transport of the spherical Ceo molecule, and three differ- 
ent short-length (5,5) carbon nanotubes Cioo, Cigo and C26O along 
the zigzag nanoribbon of the size 157 x 2 nm 2 under the action of a 
heat flow. Shown is the evolution of the coordinate of the molecule 
center for four cases (curves 1, 2, 3 and 4, respectively), (b) Temper- 
ature profile along the nanoribbon, a solid straight line corresponds 
to the temperature gradient of dT = 0.95 K/nm. Grey shaded areas 
show the edge areas of the nanoribbon subjected to two Langevin 
thermostats at T+ = 400 K and XI = 200 K, respectively, (c) 
Temporal evolution of the z projection of the center of the pentago- 
nal face relative the molecule's center of mass, corresponding to the 
curves 2 and 4 in (a). 



field E along the axis x, we create a constant external force 
F = (1, 0, 0)eE, where e is the electron charge. 

Typical dynamics of Fe~@Cgo molecule is presented in 
Fig. 6. Under the action of the electric field the compos- 
ite molecule moves along the nanoribbon towards its right 
end not detaching from it. The speed of motion grows 
monotonously before reaching the maximum velocity v m . For 
the field E = 10~ 7 V/m, the maximum velocity v m = 
390 m/s which is not sufficient to overtake the effective po- 
tential barrier created by the nanoribbon edge, so the molec- 
ular complex is reflected from the right end and remains on 
the nanoribbon, as seen from Fig. 6 (a). However, for higher 
strength of the external field E — 2 x 10~ 7 V/m the molecule 
reaches the maximum velocity v m = 540 m/s, and it may drop 
off the right end of the nanoribbon. More detailed analysis 
indicates that the sliding motion is observed only at the begin- 
ning of the motion, whereas the main motion of the Fe _ @Cgo 
molecule is its rolling along the nanoribbon in the direction of 
the applied electric field, see Fig. 6 (b). 



VI. THERMOPHORESIS OF C 60 MOLECULES 

In order to model the motion of a spherical Cgo molecule 
under the action of a constant heat flow, we consider the 
graphene nanoribbon of the size 157 x 2 nm 2 and place its 
edges into the Langevin thermostats at different temperatures: 
first 20 sections - at the temperature T + = 330 K, and the 
last 20 sections, at the temperature T_ = 270 K. After some 
transition period, we observe the formation of a constant tem- 
perature gradient dT = 0.30 K/nm, as shown in Fig. 7 (b). 
Then, we place a molecule Cgo on the 22-nd section of the 
nanoribbon and study its temperature-driven dynamics. 

Our modeling demonstrates that Cqo molecule moves to- 
wards the right end of the nanoribbon under the action of 
the temperature gradient with the constant velocity v m = 
190 m/s, see Fig. 7 (a). At the initial stage, the typical motion 
of the molecule is sliding, but after some time the molecule 
starts rolling along the nanoribbon, see Fig. 7 (c). 

For larger temperature difference, e.g. T + = 400 K and 
T_ = 200 K, the temperature gradient is larger, dT — 
0.95 K/nm [see Fig. 8 (b)], and the molecule speeds up reach- 
ing the maximum velocity v m = 440 m/s, and then it, is 
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reflected elastically from the edge of the nanoribbon, see 
Fig. 8 (a). 

VII. THERMOPHORESIS OF CARBON NANOTUBES 

To compare the transport properties of fullerenes and car- 
bon nanotubes and evaluate the effect of shape on ther- 
mophoresis, we consider the mobility of short-length (5,5) 
carbon nanotubes Cioo, Ciso and C260 These nanotubes have 
the form of short cylinders with the diameter D = 0.70 nm 
(that equals to the diameter of the Ceo molecule) but with dif- 
ferent lengths L = 1.18, 2.19 and 3.17 nm, respectively [cf. 
Fig. 1(a) and Figs, l(b-d)]. The interaction energy will be 
maximal provided the nanotube is place along the nanorib- 
bon, so that rocking is hardly possible because for this lat- 
ter case the nanotube should be placed across the nanorib- 
bon. Therefore, we anticipate that the most typical motion 
for short-length nanotubes would be their sliding along the 
surface of the nanoribbon. This is indeed confirmed by our 
numerical modeling. 

Again, we consider the graphene nanoribbon of the size 
157 x 2 nm 2 and place its edges into the Langevin thermostats 
at different temperatures: first 20 sections - at the tempera- 
ture T + = 400 K, and the last 20 sections, at the temperature 
T_ = 200 K. After some transition period, we observe the for- 
mation of a constant temperature gradient dT = 0.95 K/nm, 
as shown in Fig. 8 (b). Then, we place a nanotube next to the 
left end of the nanoribbon and study its temperature-driven 
dynamics. 

Our modeling demonstrates that the nanotube moves to- 
wards the right end of the nanoribbon under the action of 
the temperature gradient with the constant velocity v m , see 



Fig. 8 (a). The maximum velocities are: v m — 444 m/s, 
for Cqo molecule, and v m = 396, 328 and 283 m/s, for the 
nanotubes C100, Ciso, and C26O1 respectively. Moreover, our 
study indicates that the nanotubes do not roll but instead they 
slide along the surface, this is indicated by the study of the 
temporal evolution of the z projection of the center of the pen- 
tagonal face relative the nanotube's center of mass, as shown 
in Fig. 8 (c). 



VIII. CONCLUSIONS 

We have studied the motion of fullerene molecules and 
short-length carbon nanotubes along the deformable surface 
of a graphene nanoribbon in several cases, including the 
stochastic motion and directed motion in the presence of an 
electric field or a thermal gradient. We have revealed that the 
mobility of carbon structures depends substantially on their 
shape and symmetry, so that larger molecules and more asym- 
metric structures demonstrate much lower mobility than Cqo 
molecules. We have demonstrated that both electorphoresis 
and thermophoresis can be employed to control carbon struc- 
tures and fullerenes through their shape, size, and possible 
inclusions. 
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